A note on observables for counting trails and paths in graphs 
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We point out that the total number of trails and the total number of paths of given length, 
between two vertices of a simple undirected graph, are obtained as expectation values of specifically 
engineered quantum mechanical observables. Such observables are contextual with some background 
independent theories of gravity and emergent geometry. Thus, we point out yet another situation in 
which the mathematical formalism of a physical theory has some computational aspects involving 
intractable problems. 



I. INTRODUCTION 

When the length is part of the input, counting trails 
and paths in graphs is usually an expensive computa- 
tional task. For example, counting the number of Eule- 
rian trails and Hamiltonian cycles are #P problems (see 
@ and [Hj], respectively). 

In this note, we point out that the number of trails and 
the number of paths having generic given length can be 
obtained as expectation values of specifically engineered 
quantum mechanical observables. The observables arise 
from an operational construction used to associate en- 
ergy to graphs in certain background indepedent models 
of gravity (see [HQ)- The Hamiltonian of the system 
depends only on minimal information encoded in graphs, 
like, for example, the degrees sequence and the length 
of trails. The induced dynamics tends to maximize the 
number of trails of certain preferred lengths during the 
time evolution of the system. It has been argued that 
such models can exhibit a phase which describes a emer- 
gent geometry. We are interested in remarking some com- 
putational aspects of a physical model and it is out of our 
scope to propose any algorithm for ffP problems. 

The remainder of this article is structured as follows: 
in Section II we give the necessary definitions. Section III 
considers trails; Section IV, paths. Section V concludes 
the paper with a brief discussion. 



II. COUNTING TRAILS AND PATHS 

Let G = (V, E) be a simple undirected graph, where 
V(G) = {1, n}. A walk of length I in G is a non-empty 
sequence of vertices V1V2 ■ ■ ■ ViVi+i, such that {u,, G 
E(G), for every i < I. The vertices v\ and vi+\ are said 
to be the initial and final vertex of the walk, respectively. 
If Vi = vi + \ then the walk is said to be closed. If Vi ^ 
then the walk is said to be open. When this is the 
case, we chose to omit the adjective. A trail is a walk in 
which no edges are repeated. A Eulerian trail is a trail 
of including all edges. A graph with a Euler trail is said 



to be Eulerian. A path is a trail in which no vertices are 
repeated. A cycle is a closed path. A Hamiltonian cycle 
is a cycle of length n, that is, a cycle including all vertices 
of G. A graph with a Hamiltonian cycle is said to be 
Hamiltonian. See the book by Diestel for a reference 
on the concepts and terminology of graph theory. 

The adjacency matrix of G is a binary n x n matrix, 
denoted by A(G), with A(G) U<V = 1 if {u, v} G E(G) and 
A(G) U _ V = 0, otherwise. Let w(G,l,u,v) be the number 
of walks of length I in G, with initial vertex u and final 
vertex v. It is well-known that w(G,l;u,v) = A(G) UV , 
for all u, v G V(G), even if u = v. During our discussion, 
it is useful to define a formal adjacency matrix A(G), by 
replacing each A(G) u ,v — 1 with an independent variable 
e u<v , where [e UiV ,e WtZ ] — 0, for all u,v,w,z G V(G). For 
instance, for the 4-cycle C4, 



A(Ci) = 



/0110 
10 1 
10 1 
V 1 1 Q 



and 



A(C 4 ) 



( e i2 e i3 

e 2 i e 2 4 

e 3 i e 34 

V e 42 e 43 



The walks of lengths 3 between vertex 1 and vertex 2 are 
given by A{Ci)\ 2 = ef 3 ei, 2 + ei, 2 e| 4 + ei j3 e3 i 4e 2 , 4 + e\ 2 
and then w(C 4 , 3; 1,2) = 4. We denote by t(G, I; u, v) and 
p(G, I; u, v) respectively the number of trails and paths of 
length I in G, with initial vertex u and final vertex v. 



A. Trails 

Let Ti.2 = C 2 be a two-dimensional Hilbert space with 
the orthonormal basis {|0), |1)}. Let a and be opera- 
tors obeying the relation 



{a,a f } 



aa) + a) a = I. 
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These operators act on {|0), |1)} as follows: a\0) — 0|0), 
at|0) = |1), a\l) = |0) and = 0|1). The Hermitian 

combination N = a* a is the so-called particle number 
operator and it has the property that 

N\k) = o)a\k) = k\k), 

for k = 0, 1. Each Z-th (/ > 2) normally ordered power of 
N gives : N l : \k) = 0\k). For instance, 

: N 2 : \k) = aMaa\k) = 0\k). 

This is equivalent to say : N l : only has zero eigenvalues. 
Let us now define the space 

^ (g) , 

u,v£V(G) 

where dimHy 2 = 2^ 2 \ Each pair {u, v} is associated to 
a space 

( n ^)u,v= span{|0„, v ), \l u ,v)}- 

All vectors 

\ x ) = (g) \ x u,v), 
u,veV(G) 

being x u>v G {0, 1} and a; G {0, 1}( 2 ) (i.e., the set of bi- 
nary strings of length n(n — l)/2) form an orthonormal 
basis of Hy 2 - The state in Hy 2 corresponding to the 
graph G is the basis state \iPg), m which \x u-v ) = |1) if 
{u, v} G E(G) and |0), otherwise. The state \ipc) needs 
n(rt — l)/2 qubits to be encoded. Operators acting on 
the space Hv 2 can be defined by making use of the op- 
erators a and acting on the individual spaces (H.2)u,v 
In particular, it is possible to define number operators 
acting on each copy of {H%)u,v These operators are de- 
noted by N u<v and are defined by N u>v \x) = x UlV \x), for 

x G {0, 1}( 2 ). Each operator N u v returns the occupation 
number x u _ v . More explicitly, N UjV is defined by 

N UtV \x) = l(g)l(g)---(8)l(8)a t a(8)l(g)---® l|ac> = x U)V \x). 

Since the operators N UfV act on different subsystems of 
Tlv 2 , they all commute with each other. 

Now, let us define a matrix N with elements 

»j f N u>v , ifu ^ v; 
u ' v ~\0, otherwise. 

Note that the matrix N is not an operator on Jiyi in 
the usual sense; it is the elements of N that act naturally 
on 7iy2. Thus, the action of N on a state should be 
understood as occurring entry-wise. The evaluation of 
the matrix N, using the state \4>g), gives the same entries 
of the adjacency matrix, that is, 

{ipa\N u , v \ip a ) ■= (MNuMg) = A{G) U<V . 



Powers of N can be defined recursively as N = 
NN'^ 1 . The entries of N' are sums of sequences of num- 
ber operators acting on different subspaces. By consid- 
ering the expectation values of the normally ordered en- 
tries, we can state the following proposition: (Recall that 
a trail is a walk in which no edges are repeated.) 

Proposition 1 Given a graph G, the total number of 
trails of length I inG, with initial vertex u and final vertex 
v is 

t(G,l;u,v) = (iJ G \:-N l u>v :\^G). (1) 

Proof. (Sketch) Keeping in mind the formal adjacency 
matrix A{G), we need to point out the following two facts 
only: 

• If a product of number operators in one of the terms 

v contains a member acting on a pair of non- 
adjacent vertices then the operator yields zero and 
that term does not contribute to the expectation 
value. 

• When a term in v contains more than one copy 
of a number operator acting on a particular edge, 
that term also does not contribute to the expecta- 
tion value, because : N l : \k) — 0\k). 

It follows that the only combinations of number op- 
erators that do not give a vanishing contribution to 
t(G, l,u,v) correspond to trails with initial vertex u and 
final vertex v. Thus, Eq. ([TJ) gives this number of trails as 
an expectation value of a quantum mechanical operator. 



Notice that : N™ u : counts the number of Eulerian 
trails in G, if \E(G)\ = m. 

It is noteworthy to remark that the logarithm of the 
dimension of the space Tiy 2 is polynomial in the number 
of vertices. Also, notice that the operators v are in- 
dependent of G. For the purpose of counting trails in a 
specific graph, a similar procedure may be applied taking 
a Hilbert space 

He ■= (g) (W2)«,», 

{u,v}eE(G) 

where dimTi.E = 2\ E ( G ^. In this case, the matrix of 
operators N needs to be slightly modified so that some 
of its entries are zero from the beginning, rather than 
number operators acting on empty states. 

Proposition [1] is based on the equation : N l : \k) = 
0\k). It may be worth remarking that there is way of 
counting trails without making use of normal ordering. 
This can be done by defining the matrices 

n _f a u>v , ifu^v; 
u ' v ~ \ 0, otherwise. 
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It is in the same spirit of matrix N Mjt) , but the entries 
are annihilation operators, a u v , rather than number op- 
erators N UfV . This matrix of operators is not Hermitian, 
i.e., ^ D. Here the Hermitian conjugate is defined 
entry-wise by (D^) u ,v = (Du,v) ■ It can be shown that 

t(G,l;u,v) = (ip G \:N l UiV : \$ G ) 

B. Example 

We consider the 4-cycle C 4 : the set of vertices is 
V{C 4 ) = {1,2,3,4}; the set of edges is E(C A ) = 
{{1,2}, {1,3}, {2, 4}, {3, 4}}. Since \V(C 4 )\ = 4, we asso- 
ciated to this graph an Hilbert space Hv 2 of dimension 
2 6 = 64. Specifically, 

H V 2 = {H 2 ) h2 ® (W 2 )i, 3 ® ■ ■ ■ ® (W 2 ) 3; 4 • 

The state associated to C 4 is 

|Vg) - 1110011). 

Regarding the operators iV U)t ,, we have, for instance 

7Vi j2 |110011) = 1|110011) 

and 

7Vi j4 |110011) = 0|110011). 

In fact, 

(110011|TVi, 2 |110011) = A(C 4 )i, 2 = 1 

and 

(110011|JVi,4|110011> = A(C 4 ) 2 ,4 = 

Observe that A{Ci)\ 2 — e? i3 ei i2 + ei i2 e2 ;4 + ei j 3e3 ; 4e 2i 4 + 
e\ 2 - For the vertices 1 and 2, we have, 

(Vg|: N? >2 : |Vg) = (V> G |: ^3^,3: |Vg) 
+ (Vg I : N 1>2 Ni A : |Vg> 

+ (^ch ^1,3^3,4^2,4: \i>c) 

+ {4> G \: Nf X . |Vg) 

= (V>g|: ^1,3^3,4^2,4: IV-g) 
= 1 

Indeed, £(C 4 , 3; 1, 2) = p(C 4 , 3; 1, 2) = 1, something that 
: Nf 2 ■ is able to detect even if u>(C 4 , 3; 1, 2) = 4. 

C. Paths 

In this section, our working space is TLv '■— 
®»eF(G)(^2)»i where dimTLv = 2™, given that 



|y(G)| = n. This is the space usually assigned to net- 
works of spin 1/2 quantum mechanical particles. The 
space (Ti,2) v is associated to the vertex v and (TL 2 ) V = 
span^), |1„)}. All vectors \x) = ® veV { G ) \x v ), being 
x v G {0, 1} and x £ {0, 1}", form an orthonormal basis 
of TLv- The state in TLv corresponding to the graph G 
is the basis state |ll...l) = n |1). As we have done in 
the previous subsection, we can define number operators 
acting on each (TL 2 )v These operators are denoted by 
N v . 

Let M be a matrix of operators with entries defined as 
follows: 

!4ax(«,»), if{u,w} £ E{G) and v > u; 
N min{u , v ) , if {u, v} e E(G) and u > v; 
0, otherwise. 

As in the previous subsection, the elements of M are still 
number operators, but this time acting on TLv- However, 
unlike the adjacency matrix, M is not symmetric. Powers 
of M arc defined via matrix multiplication, for example, 
by the recursion M' = MM' -1 . Recall that a path is a 
trail in which no vertices are repeated. 

Proposition 2 Given a graph G, the total number of 
paths of length I in G, with initial vertex u and final 
vertex v is 

p(G,l;u,v) = (ll...l\:M l UiV : |11...1>. 

Proof. (Sketch) The operators Mj, „ are sums of prod- 
ucts of number operators. Since the entries in M arc 
nonzero only in the positions where the adjacency ma- 
trix is nonzero, each of the terms in M^ „ also corre- 
spond to walks. Due to the normal ordering conven- 
tion, a term in which the same vertex appears more than 
once does not contribute to the expectation value. Thus, 
(ip G \ : M. l u v : \<p G ) has the desired interpretation. ■ 

A similar construction could be made up with a "sym- 
metric version" of M. 

With respect to the Hilbert space TLv, we define a 
matrix of operators F with entries 

{a ma x(«,«) if {u, v} e E(G) and v > u; 
a m in(u,v) if { u , v} e E(G) and u > v; 
otherwise. 

Notice that F u ^ v depends on E(G) and thus on the graph 
G. Then, the transition amplitude 

p(G,Z;u,u) = <00...0|Fi Mi |ll...l> ) 

is the number of cycles of length I containing the vertex 
u. Taking I = n, we have p(G, n; u,u) > if and only if 
the graph G used to construct F Ui „ is Hamiltonian. 
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III. CONCLUSION 

For the purpose of counting trails and paths, we have 
described quantum observables in Hilbert spaces whose 
logarithm of the dimension is a polynomial in the num- 
ber of vertices. The numbers of trails and paths can 
be obtained as expectation values of these observables. 
The states involved are the pure states |ll...l) and \ipc)- 
These states can be prepared efficiently. While it is clear 
that the observables raise no issues about uncertainty, it 
is most likely that these are not efficiently implementable 
in a quantum computer (e.g., by phase estimation). The 
reason behind this thought is based on the fact that our 
observables are used to approach #P problems. Despite 
this, it may still be instructive to describe their form for 
special classes of graphs, to determine their complexity, 
and to describe the physics required for the implementa- 
tion. 

It is valuable to point out that the literature contains 
so far a number of examples of quantum observables for 
solving computational task, which are either not known 
to be efficiently implementable or implementable with 
poor accuracy. Among these, observables for t he g raph 
isomorphism problem [5[ and for the permanent [111 ] (see 
also [13|). It has also be shown that some mathematically 
well- defined observables allow to solve the halting prob- 
lem, in contradiction with the Church- Turing thesis Q. 
Recall that the Church- Turing thesis asserts that every 
function which can be computed by what we would nat- 
urally regard as an algorithm is a computable function, 
and viz. 



In the present notes, we have highlighted yet another 
situation in which the mathematical formalism of a phys- 
ical theory has some computational aspects involving in- 
tractable problems. On the other hand, the computa- 
tional complexity side may suggest some potential limita- 
tions on the physical picture (see [l| , for a detailed survey 
on these ideas) - that is, the physical picture is simply not 
plausible. - With the same perspective of [8j, it is legiti- 
mate to study the properties of background independent 
models of gravity and emergent geometry as computa- 
tional devices (as cellular automata, for instance). In- 
deed, in the model studied in Q , the Hamiltonian of the 
system keeps track of trails of given lengths. Still this is 
inline with natural phenomena of self-organization at the 
microscopic scale. For example, the folded 3-dimensional 
conformation of a protein is believed to be its lowest free 
energy state. The 3-dimensional models describing the 
folding process of a protein as free energy minimization 
problems are TVP-hard (see, e.g., Q). Finally, since the 
entries of the observables described in this paper can 
be seen as Grassmann numbers (because corresponding 
to fermionic operators), the observables may have some 
analogy to certain matrix model (see, e.g., [Io|'). where 
the partition function is given by their weighted sums. 
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